Molecular production in two-component atomic Fermi gases 
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We provide a practical approach to the molecular production via linear downward sweeps of Feshbach res- 
onances in degenerate Fermi gases containing incoherent mixtures of two atomic spin states. We show that 
the efficiency of the association of atoms is determined just by the Landau-Zener parameter in addition to the 
density of the gas. Our approach of pairwise summation of the microscopic binary transition probabilities leads 
to an intuitive explanation for the observed saturation of the molecular production and recovers all atomic loss 
curves of C.A. Regal et al. [Nature (London) 427, 47 (2003)] as well as K.E. Strecker et al. [Phys. Rev. Lett. 91, 
080406 (2003)] without adjustable parameters. 
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Ever since the achievement of Bose-Einstein condensation 
in a dilute vapour of ^^Rb atoms, there has been considerable 
interest in extending the range of species for studies of cold 
gases. Several experiments | lll2,bl|4|,|3lal2llsll2ll have now 
demonstrated the production of cold molecules via magnetic 
field tunable Feshbach resonances. One of the techniques re- 
lies upon the adiabatic association of pairs of cold atoms to 
molecules through a linear sweep of a Feshbach resonance 
level across the zero energy threshold of the colliding atoms. 
This technique has been applied to atomic Bose-Einstein con- 
densates |4, 5, 9], as well as degenerate two-component Fermi 
gases 1 3, 6] containing a balanced incoherent mixture of Zee- 
man levels of a single atomic species. In contrast to cold 
Bosons, degenerate Fermions inevitably occupy a range of 
single particle energy levels up to the Fermi edge. There 
have been concerns that the spread in collision momenta of 
Fermionic pairs could prevent an efficient molecular produc- 
tion in an adiabatic sweep. Despite this intrinsic difficulty, a 
significant molecular conversion of up to approximately 50% 
of the atoms was reported in Refs. |3, 6], and higher conver- 
sion efficiencies have been reported when much slower ramp 
speeds are used Hli^l. The few attempts O El El Q 
to explain these observations are all based on rather involved 
field theoretic models. While Refs. |12, 13] provide simula- 
tion methods for the molecular production in the presence of 
Fermi seas, Ref. 1 11] suggests the maximum molecular pro- 
duction efficiency of 50% in Ref. 1 3J to be a consequence of a 
significant temperature dependence of the conversion process. 

In this letter we derive a practical universal formula for the 
molecular production via a linear sweep of the magnetic field 
strength B in a degenerate two-component Fermi gas in the 
short time limit where only a single collision of a given atom 
occurs during the sweep. Our analytic treatment relies upon a 
pairwise summation of the probabilities for the atomic associ- 
ation, which naturally explains the fast ramp observations of 
Refs. l2iy]- We show that the molecular production is largely 
insensitive to the cold collision momenta and depends only 



on the product ObgAB of the background scattering length Obg 
and the resonance width AB, in addition to the atomic mass 
m, the ramp speed B as well as the density of the gas. Despite 
its simplicity, our analytic formula leads to predictions for the 
molecular production that are virtually indistinguishable from 
the observations of Refs. |3, 6] within the experimental error 
bars. 

We shall first discuss the association of atomic pairs in the 
absence of the surrounding gas. A homogeneous magnetic 
field B splits the atomic hyper-fine angular momentum states 
into a set of Zeeman levels (f,mf), labelled by the projec- 
tion quantum number mj of the total atomic spin along the 
axis of the magnetic field and the / value with which it cor- 
relates adiabatically at zero field. Indistinguishable Fermions 
with equal spins are not subject to s-wave scattering and do 
therefore not contribute to molecular formation. We shall thus 
consider pairs of atoms whose internal states correspond to 
different Zeeman levels at asymptotically large inter-atomic 
separations. We denote the .s-wave binary scattering channel 
of such a pair of unlike Fermions as the entrance channel. 

The physical concept of associating a pair of cold atoms is 
closely related to the variation of its energy spectrum under 
adiabatic changes of B. Figure ^ shows a coupled channels 
calculation of these energies for the example of the interact- 
ing pairs of "*"K atoms in the (/ - 9/2, ot/ = -9/2) and 
(/ = 9/2, m/ = -5/2) Zeeman states of Ref. |3]. The fig- 
ure shows that when B varies from high to low fields across 
the zero energy resonance at Bo = 224 G (1 G= 10""* T), a 
molecular bound state (pbiB) appears whose energy E\,{B) is 
indicated by the solid curve. The zero of the binding energy 
E\,(B) coincides with the magnetic field strength Bq at which 
the scattering length has a singularity. 

The strong magnetic field dependence of the cold inter- 
atomic collisions is due to the near degeneracy of the energy 
EresiB) (dotted dashed line in Fig.Q of a single closed channel 
vibrational state (/i,es (the Feshbach resonance level) with the 
threshold for dissociation in the entrance channel (dotted line 
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FIG. 1: The binding energy E\,(B) (solid curve) of the highest ex- 
cited dressed vibrational bound state of a pair of """K atoms versus 
the magnetic field strength B in the vicinity of the 224 G zero en- 
ergy resonance |3]. The zero of energy (dotted line) is chosen, for 
each magnetic field strength, as the Zeeman energy of a pair of dis- 
tinguishable Fermions at asymptotically large spatial separation. The 
zero of the dressed state energy E\, coincides with the singularity of 
the scattering length in the entrance channel at the magnetic field 
strength Bg. The dotted dashed line indicates the energy E^^siS) of 
the bare Feshbach resonance state. The dashed line corresponds to 
the energy E_i of the highest excited bare vibrational level of the 
entrance channel potential. 



in Fig.0. This allows us to use a two-channel configuration 
interaction representation in terms of the bare basis 

states, i.e. the states in the absence of inter-channel coupling. 
These bare states include the Feshbach resonance level (pies 



and the entrance channel continuum states 



associated with 



the collision momentum p. The kinetic energy p^/m as well as 
the magnetic moment fi^es - dEres/dB of the resonance level 
are independent of B. Consequently, the determination of the 
probability for the association of an atom pair with the relative 
momentum p can be described by the Landau-Zener approach 
lUsillill . We note, however, that although we have chosen a 
bare state representation, the bound molecules produced in a 
linear sweep of the magnetic field strength are described by 
the dressed state (f>b{B), which accounts for the inter-channel 
coupling. 

The relative momenta in a zero temperature two-component 
Fermi gas are limited by the maximum of the Fermi momenta 
of each component, whose associated energy scales are typi- 
cally on the order of fiK (in units of the Boltzmann constant). 
These energies are too small to be resolved in Fig.^ In the 
limited range of cold collision momenta the probability P for 
the association of a single pair of Fermions in a large sample 
volume is isotropic, i.e. it depends just on the wave number 
k = p/h associated with the collision momentum p. Neglect- 
ing phenomena related to the presence of the surrounding gas, 
considerations described in Ref. IJdl lead to the formula L17.1 : 
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The inter-atomic interactions enter the probability in terms of 
the Landau-Zener parameter 1 15, 16]: 
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Here W{r) is the potential describing the inter-channel cou- 
pling. We note that {(^reslW'l^p^') and consequently 5lz are 
independent of k provided that |flbgl is sufficiently small to sat- 
isfy the relation k\a\,g\ <k 1. The k dependence of Eq. Q is 
due to downward transitions into lower energetic continuum 
states that compete with the molecular formation in an iso- 
lated two-body system Lisll . 

We shall show in the following that the molecular produc- 
tion in the limit of fast magnetic field ramps as well as in the 
opposite saturation limit can be determined in a straightfor- 
ward manner, provided that multiple collisions of a given atom 
can be neglected. To this end, we shall first consider the limit 
of ramp speeds sufficiently high that Eq. Q is well approx- 
imated by first order perturbation theory in the inter-channel 
coupling. A linear expansion of Eq. Q in 6lz then shows 
that P(k) is independent of k, and reduces to the linearised 
Landau-Zener formula for a single curve crossing 



P = 2n6Lz- 



(3) 



Despite the fact that the association of a given atom pair in 
the large sample volume 'V is a rare event, the molecular pro- 
duction in an incoherent mixture of Fermions can be rather 
efficient because each of the A^i atoms in one of the spin states 
has all N2 atoms of the other component to interact with. In 
accordance with classical probability theory, the total num- 
ber A^b of diatomic molecules in the state 0b is given by the 
pairwise incoherent sum of the microscopic transition proba- 
bilities: 



A^b = N1N2P 



(4) 



Here A^iA^2 is the number of pairs of unlike atoms. Equations 
(|3} and then lead to the following fraction of atoms asso- 
ciated to molecules: 

2A^b /N = 2a'( 1 - a)NP = 2q'( 1 - a) [Inn'VdLz] ■ (5) 

Here n = N/'V is the local density and a is the minimum of 
the fractions A^i /N and N2/N. Since 'Vc^lz is independent of 
'V, Eq. (|5} depends just on the local density of the gas. The 
pre-factor a(l - a) accounts for the size of the sub-ensemble 
of pairs of unlike Fermions and the factor of two indicates the 
number of atoms per molecule, while Imi'VSuz can be inter- 
preted as a density dependent Landau-Zener transition proba- 
bility OEl- 

The pairwise summation of microscopic transition proba- 
bilities is justified when the number of molecules produced 
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is small in comparison to the total number of atoms 
In this limit of high ramp speeds the molecular production 
does not crucially depend on the presence of the Fermi seas. 
As the ramp speed decreases, Eq. Q becomes momentum 
dependent, since, unlike for fast ramps, there can be decay 
during the ramp to lower energy continuum states of the two- 
body system with £ > 0. In the presence of the surrounding 
two-component Fermi gas, however, these downward transi- 
tions would lead into occupied modes and are therefore pre- 
vented. In order to extend Eq. (|5} to lower ramp speeds, 
we shall thus follow the approach of Refs. Iisl \l6k and re- 
place the density dependent linearised Landau-Zener formula 
Itth^Slz by the complete Landau-Zener transition probabil- 
ity 1 - exp(-27rn^6i,z) for a linear curve crossing in the ab- 
sence of downward transitions to other continuum states. The 
corrections to this procedure involve the dynamic depletion 
of the available pairs during the atomic association. Com- 
parisons between genuinely many-body approaches and the 
Landau-Zener formula 1 16| for the case of a condensate of 
identical Bosons show that the dynamic depletion of the pairs 
changes the functional dependence of the molecular produc- 
tion efficiency on the ramp speed. The absolute magnitude of 
these corrections, however, is typically small in comparison 
with the total depletion of the atoms and can be neglected. The 
molecular production efficiency in a two-component Fermi 
gas may thus be approximated by: 



2Nb 
N 



= 2a(l-a) 



1 - exp -n 
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The molecular production efficiency can not exceed 50%, 
given the ansatz in Eq. This limit is quite intuitive as, 
given only one chance to associate, a Fermion in the larger one 
of the two components can find an interacting partner with a 
probability of 1/2 at best. The saturated efficiency in mixtures 
with unbalanced populations of the spin components can only 
be lower and is given by 2a(l - a). The exponential func- 
tion in Eq. (|6} provides a smooth interpolation between the 
straightforward asymptotic limits of the ramp speed. 

The range of validity of our approach is set by the assump- 
tion that only a single collision of a given atom can occur 
during the association process. When the ramp speed is cho- 
sen sufficiently low for multiple collisions to become signif- 
icant, efficiencies exceeding 50% can result |7, 10| from the 
same physical mechanisms that can convert an atomic gas to a 
molecular gas at constant magnetic field strength |8, 19, 20]. 
To this end, the experiments of Ref. |10] applied magnetic 
field ramps at least ten times slower than the slowest ramp 
of Ref. |3] for Feshbach resonance parameters and densities 
comparable to those of Ref. 1 3 1 . 

To compare the predictions of Eq. (|6j to the experimental 
results of Refs. IHl^, we have applied the local density ap- 
proximation, i.e. we have averaged Eq. (|6j in accordance with 
the density distribution of the trapped gases. We have deter- 
mined the density profiles from the zero temperature Thomas- 
Fermi approximation |21], adapted to the experimental trap 
geometries. Figure 121 shows the number - 2N^ of free "^''K 
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FIG. 2: The number of atoms remaining versus the inverse ramp 
speed 1/5 for a linear sweep of the magnetic field strength across 
the 224 G zero energy resonance in a two-component Fermi gas of 
'*''K. The experimental data of Ref. 1 3] (circles) are compared to the 
results of the pairwise summation approach of Eq. ^6} with a = 1/2. 
The solid and dashed curves indicate local and mean density approx- 
imations, respectively. 



atoms remaining in the experiments of Ref. f3] as compared to 
the predictions of Eq. (|6|l. We have estimated the total number 
of atoms to be = 1.25 x 10^ (cf. Fig.|2j and used the back- 
ground scattering length Obg - 174 aq 1 221 (flo = 0.052918 
nm), the resonance width AB = 9.7 G |22], and the radial 
angular trap frequency Wr = 2?! x 215 s"', which, together 
with the aspect ratio of 70 L2], largely recover the experimen- 
tal peak density of Wpk = 1.4 X lO'^ cm"^ (3j]. This procedure 
yields a mean density of « = 9 x 10'^ cm"^. The predictions 
of the mean and local density approaches to Eq. clearly 
follow the slope as well as the magnitude of the experimental 
data. We obtain a similar agreement with unpubhshed lower 
density (wpk = 9 x 10^^ cm"^) data of Regal et al. 

To further confirm the validity of Eq. (|6}, Fig. |3 shows 
an analogous comparison with the results of Ref. 16]. In 
these experiments the molecules were produced by a linear 
sweep of the magnetic field strength across the 543 G zero 
energy resonance of ''Li in the entrance channel associated 
with the set of quantum numbers (/ - 1/2, my = +1/2) 
and (/ = 1/2, TO/ = -1/2) of the interacting atoms. We 
have performed an exact coupled channels calculation to de- 
termine the resonance parameters. From these considerations 
we obtain |flbgAB| - 5.916 aoG. The Thomas-Fermi density 
profile used in the local density approach is determined by 
the Fermi energy E^, which in the case of the combined har- 
monic and box-Uke trapping potential of Ref. is given by 

= [l5nNn^cx?J{i V2^L)]^^^ Here w,- = 27rx800 s"' is the 
angular frequency of the radial harmonic trap and L - 480 /vm 
is the size of the axial box-like potential 1^. This leads to a 
mean density of n = 4 x lO'^ cm"^. 

The agreement with the experimental data in Figs. |2] and 
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FIG. 3: The fraction 1 - 2Nb/N of remnant atoms versus the inverse 
ramp speed 1/5 for a linear sweep of the magnetic field strength 
across the 543 G zero energy resonance in a two-component Fermi 
gas of *Li. The circles indicate the experimental data of Ref. |6]. The 
solid and dashed curves are predictions of Eq. J6j with a = 1/2 in 
the local and mean density approximations, respectively. 



13 suggests that the pairwise summation approach recovers all 
the relevant aspects of the physics involved. This implies that 
the molecular production efficiency is determined by the den- 
sity of the gas, the ramp speed and the atomic mass, in addi- 
tion to a single universal parameter of the resonance enhanced 
colhsion physics, the product of the background scattering 
length and the resonance width. The density profile provides 
a direct dependence of the molecular production on the Fermi 
statistics. Our results indicate, however, that even a mean 
density approximation yields a rather useful estimate of the 
number of molecules produced. Our ansatz provides an ex- 
cellent approximation for the molecular production efficiency 
in the limit of very fast downward sweeps of the Feshbach res- 
onance level in two-component Fermi gases and gives a 50% 
efficiency as a natural limit, provided that phenomena related 
to multiple collisions of a given atom are negligible. The pair- 
wise summation approach has a rather wide range of applica- 
bility in dilute gases, including cold Bosons ITa flH fisll and 
Fermions as well as their mixtures. The main requirement is 
the diluteness of the gas in combination with a binary transi- 
tion, which is sufficiently fast that multiple collisions can be 



neglected. Such a binary transition may be realised by a single 
linear sweep of the magnetic field strength, but it can as well 
involve fast sequences of magnetic field pulses [18] applied in 
the experiments of Refs. ll),y]. 
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